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14.  Abstract 


This  paper  was  presented  during  the  42nd  Meeting  of  the  Structures  and  Materials  Panel  in 
Ottawa  in  April  1976. 

It  deals  with  a serious  difficulty  in  unsteady  aerodynamics,  that  is  the  prediction  of  th.j 
pressure  field  induced  by  the  rotation  of  a control  surface.  Much  work  has  already  been 
done  on  this  subject  in  subsonic  flow,  but  this  is  one  of  the  first  approaches  to  the  super- 
sonic problem.  Predictions  have  been  made  by  two  methods  developed  separately  by  BAC 
and  MBB.  They  have  been  compared  with  windtunnel  tests  made  at  NLR  us»ng  more  than 
80  pressure  tubes.  Pressure  distributions,  hinge  moments  and  lift  have  been  measured  for 
different  sections  of  the  wing.  As  the  two  theories  that  have  been  used  arc  linearised,  the 
agreement  between  theory  and  experiments  is  not  perfect  but  appears  to  be  adequate  for 
flutter  speed  prediction. 
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This  paper  by  Lodge  and  Schmid  was  presented  to  the  Sub-Committee  on  Aero-  i 

elasticity  and  Unsteady  Aerodynamics  cn  5 April  1976  during  the  4,?nd  Meeting  of  the  J 

Structures  and  Materials  Panel  in  Ottawa. 

| 

It  deals  with  one  of  the  m^st  serious  difficulties  in  unsteady  aerodynamics,  that  is  ' 

the  prediction  of  the  pressure  field  induced  b>  the  rotation  of  a control  surface.  Much 
work  has  already  been  done  on  this  subject  in  si.bwnic  flow  but  this  publication  presents 

one  of  the  first  approaches  to  the  supersonic  problem.  t 

i 

Predictions  have  been  made  by  two  distinct  me'hods  developed  separately  by  BAC  ; 

and  MBB.  They  have  been  compared  with  windtunnel  tests  made  at  NLR  using  more 
than  80  pressure  tubes.  Pressure  distributions,  hinge  moments  and  lift  have  been  measured 
for  different  sections  of  the  wing. 

As  the  'wo  theories  that  have  been  used  are  linearised,  the  agreemert  between  theory 
and  experiment  is  not  perfect  but  appears  to  be  adequate  for  flutter  speed  prediction. 

This  paper  is  of  great  interest  to  the  acroelasticians  of  the  NATO  community  and 
will  help  both  flutter  prediction  and  active  control  design.  | 
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Suaaarr : 

Most  aircraft  flutter  probleas  have  featured  control  surfaces,  and  It  Is  necessary  that 
unsteady  aerodynaalc  forces  generated  by  their  notions  should  be  accurately  predicted. 
Therefore,  theoretical  and  experlaental  studies  on  a planfora  with  part-apan  control 
surface  oscillating  In  the  control  surface  rotation  node  at  loe  supersonic  Mach  nuabera 
have  been  effected  ant  the  results  of  these  are  presented  and  dlacuased.  It  la  shown 
that  these  studies  aust  be  of  a high  accuracy  so  that  the  aore  critical  aerodynaalc 
coefficients,  such  as  binge  aoaeat  daaping,  night  be  deteralned  with  confidence. 


1 . ROTATION 
Syabols 

5(F)  local  chord  of  aerofoil  (Including  control  surface) 

c aean  chord 

f frequency 

M Mach  nuaber 

p pressure 

s seal - span 

V free-stresa  velocity 

e non-dlaenslonal  jpwaer  amplitude 

s,y,s  non-dlaenslonal  rectangular  co-ord' nates  referred  to  s 

e density  of  air  * 

t n,  ,-dlaenslonal  aaplttude  of  velocity  potentiol 


Subscripts 

h . 1 . hinge  line 

1 . e . leading  edge 

t. e.  trailing  edge 

r.  rudder 

l.a.  loser  surface 

u. a.  upptr  surface 

Definitions 
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ipy" 


induced  frequency  paraaeter 
stagnation  pressure 
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AC 


- S7  */  hypvrbollc  distance 

M*-l 

^ 1 . s . **  u . a . 


A /q  * AC'  ' - 1 AC  " 

P ~ P P 


2 r f 


'"Irc.t.ar  frequency 


mzm , mm  mms 


t . m. 

l(r)  - » / fl  <1*  “ <l„  <=(y)  .Cj 

*h.  1 . 


sectlcihsl  lift 
(positive  upwards) 


My) 


t ■}. 


• 2 / ap  <x*xh.1.>  d*  ' *-  c<y> 2 *''u 


b.  1 . 


sectional  hinge  aovent 
(positive  nose  downwards) 


s t . e . _ 

s’  j j a dx  dy  * c.s.C^ 


0 x 


1 .e. 


total  lift 
(positive  upwards) 


>’  / / * p <*-*t  !,>  ix  “y  1 1. 

0 3 


• • • C„ 


*h.l. 


total  hints  aoatnt 
(positive  nose  downwards) 


2.  INTRODUCTION 


The  accurst;  of  the  control  surface  unsteady  aerodynamics  Is  ;-'rt 1 cul ar 1 y Important 
on  modern  combat  aircraft  at  low  supersonic  Mach  nuaber/high  frequency  parameter 
combinations . 

*n  order  to  assess  the  accuracy  of  current  theories  for  control  surface  unsteady 
aerodynamics,  a nomlnell)  rigid  model  has  been  designed,  built  and  tested  st  NLR. 
Amsterdam,  at  Mach  numbers  up  to  1.3  and  reduced  frequency  parameters,  based  on  se  . - 
span,  up  to  1.6.  Unsteady  pressure  distributions  Induced  over  the  main  and  control 
surfaces  by  rotntlonal  control  surface  oscillations  haco  been  measured  and  compared  with 
predictions.  Since  dynr».  tally  scaled  con* rol  surface  models  can  involve  a combination 
of  rigid  rotation  and  torilon,  which  both  ,enera*.e  basically  similar  unsteady  aerodynamic 
eifects,  model  data  fur  a i'-tatlon  mode  car  b*  used  to  quantify  tuc  accuracy  of  flutter 
calculations  Involving  the  :-eal  aircraft  vitratlon  modes. 

Predicted  a,>d  measured  results  are  pteret.ted  here  for  M * 1.1.  1.3,  and  values  of 
k from  0.7  up  to  1.1. 

3.  APPEL  DETAILS 

The  model  under  consideration  Is  a swept  main  surface  with  part-span  control  surface, 
of  aspect  rrtlo  2.0,  as  described  n ligure  la.  lit*  profile  Is  symmetric  with  a maximum 
tUl*  kness  to  chord  ratio  of  5.5%.  npt>  between  .he  main  and  control  surfaces  have  been 
kept  down  to  the  order  of  (S.l  am. 

In  the  model  construction,  both  m and  control  surfaces  were  made  as  stiff  as 
possible,  to  minimise  any  vibration  c.-/  effects.  The  main  surface  is  of  steel  snd  the 
control  wui fac»-  of  D ral , tc  reduce  the  Inertia  forces. 

Results  were  obtained  from  measurements  at  pressure  boles  e.tia.ed  along  3 atrei'-- 
wise  stations  (Figure  Is)  These  are  tore  closely  spaced  nenr  »ne  h *.'.■; e lln  and  the 
control  surface  tip. 

During  the  model  design  stagv,  research  work  Scferences  1 and  2)  revealed 
Inadequacies  In  the  calibration  factors  hitherto  applied  by  KLR  to  the  measured 
pressures,  to  account  for  tube  system  dynamics  vltn  wlnd-on.  Therefore,  special 
measures  were  token  during  these  model  tests,  to  establish  the  correct  cs 1 1 bra t ion . 
by  Installing  sooe  direct  measuring  pressure  transduc er>  (Reference  3). 

In  addition,  accelerometers  were  Installed  to  determine  the  true  model  vibration 
mode  during  osclllstory  control  surface  cxcltstion. 

4 - MNP  TUNNEL  TESTS 

Under  contract  with  B.A.C.,  Vartoti,  oscillatory  pressur.  measurements  sere  made  on 
the  model  described  above  in  the  High  Speed  Tunnel  at  the  XL*!,  Amsterdam  In  this 
closed  circuit  tunnnl,  with  s rectangular  test-section  of  2.00  metres  wide  by 
1.60  metres  high,  the  values  of  the  Reynolds  number,  based  on  mean  chord  varied  within 
the  range  of  3.5  a 10*  and  4.7  x lO*  (Figure  lb). 

The  model  was  mounted  into  the  tunnel  wnll.  sith  the  main  surface  rigidly  clamped. 

The  control  surface  excitation  equipment  waa  outside  the  tunnel  wall  Further  teats 
exhibited  that  the  repeatability  of  the  results  was  within  10%.  The  boundary  layer 
transition  was  fixed  by  Installing  strips  of  carborundum  grains  o/  nominal  size  74.  on 
both  sides  of  the  model.  These  strips  were  placed  forward  of  the  region  of  shock  onset 
(Reference  4). 


v.  - . . 





AERODYNAMIC  THEORIES 


5.1  Ou  1 1 1 r a » o t B.A.C.  Mot  hod 


The  supersonic  lifting  surface  theory  of  Sadler  and  Allen  (Reference  5)  is  a 
linearised  theory  solving  the  doublet  based  integral  equation.  This  Involves  the 
solving  of  an  Integral  equation  which  relates  the  down-wash  at  any  given  point  to 
the  p 'i turbat ion  velocity  potential. 

Integration  has  to  be  performed  over  that  portion  of  the  sing  and  the  wing  wake 
lying  in  the  Mach  fore-cone  of  that  point.  As  the  potential  la  xero  forward  of  the 
plaxfore  leading  edge  and  outside  the  wing  wako.  the  integrations  require  less 
storage  and  are  speedier  than  those  of  the  sort  popular  integrated  downwash  setbods . 

The  basic  integral  equation  !»:- 


*(*.y) 


if 

Mach 

fore-con# 


S(C,n)  K ( x - f.  , y - n > 


K(x.y)  * 


rMk  R^  1 


Onaractcr : it tc  co-ordinate*  1 and  i,  with  their  origin  at  the  pivotal  point 
(the  point  at  -thich  the  potential  is  being  calculated)  and  their  axn  parallel  to 
the  t»c  forward  pointing  Mach  Mnei,  are  lntro^i-ced.  the  linea  of  Integer  i and  3 
forming  a characteristic  a <etU  (figure  2). 

By  splitting  the  above  equation  into  its  real  and  imaginary  parts,  and  assuming 
certain  variations  of  the  real  and  imaginary  parts  of  $ along  lines  of  constant  £ 
and  known  functions  of  the  upwasft  w at  the  pivotal  point  are  found  as  weighted 
sums  of  the  potentials  at  the  mesh  points,  from  rhlch  the  potential  at  the  pivotal 
point  can  ire  determined. 

When  the  trailing  edge  is  subsonic,  the  potential  over  the  wake  region  is  needed 
and  this  may  be  calculated  from  the  condition  that  no  load  can  be  sustained  across 
,ve  wake.  Thus  we  find.- 

5 ■ : exp(-ik(x  . - x.  >)  (3) 

wake  t . e . wake  t . e . 

Where  a box  is  cut  by  the  leading  edge,  special  weighting  has  to  be  Introduced 
for  that  particular  box.  to  prevent  auy  leading  edge  inaccuracies. 

The  integral  equation  ia  »>lved  in  s "marching"  technique  such  that  every 
potential  in  the  pivotal  joints  forward  pointing  Mach  cone  is  known.  Thus,  the 


potentials  at  every  mesh  point  can  be  calculated 
matrices  are  then  derived  from  these  potentials. 


and  the  aerodynamic  work 


As  programmed  this  theory  cannot  deal  with  the  singularity  associated  with  the 
control  surface  hinge  line  incidence  discontinuity.  Tnls  problem  has  been  overcome 
b>  treating  the  whole  planform  as  a single  surface  and  smoothing  out  the  hinge  line 
discontinuity  by  a cubic  interpolation. 


5 2 Ou t 1 1 res  of  M , B B Theory 

The  unsteady  program  used  by  M B.B.  is  an  extension  of  the  characteristic  box 
method  first  proposed  by  V J E Stark  (Reference  6)  The  assumed  small  perturbation 
of  the  flow  field  implies  a lull  linearisation  of  the  potential  equation  The 
procedure  was  amended  and  programmed  by  C Rohm  (Reference  7)  and  H Schmid 
(Reference  8)  Primarily,  the  velocity  potential  of  the  harmonically  var>lng  flow 
is  evaluated  by  means  of  a source  distribution  The  basic  equation  reads  as  follows 

i « H w(x’.y')  Kfx-x'.y-y'.M.k)  dx’dy*  (4) 

Mach 

fore - cone 

The  equation  relates  the  velocity  potential  to  the  normal  velocity  w.  K 
represents  the  contribution  to  the  potential  of  a harmonically  pulsating  unit 
source  and  is  defined  in  Eq  (2) 


The  integration  region  is  limited  by  the  Mach  fore-cone  emanating  from  the 
collocation  > >lnt  and  by  the  Mach  rays,  emanating  from  the  wing  or  control  surface 
apex  Withit  this  region,  w • 0 generally  (figure  1) 


4 


Thw  turfict  in'wgral  taken  over  Mach  fore-cone  le  subdivided  by  the  program 
Into  a number  of  surface  integrals  over  rhombic  or  characteristic  boxes,  the  edges 
of  which  are  parallel  to  the  Mach  lines  of  the  fore  - and  aft  cones,  respectively. 
Within  each  box,  the  downwash  distribution  « is  assumed  to  be  constant. 

If  turn  planforn  has  a subaoolc  leading  or  trailing  edge,  or  streamvlee  parallel 
side  edges,  then  the  concept  of  diaphragms  according  to  Evvard  Is  introduced  into 
the  H.B.B.  procedure. 

Whereas  the  source  distribution  or  normal  velocity  distribution  on  the  lifting 
surface  is  known,  the  source  dlsvrlbut ion  placed  on  the  diaphragm  is  unknown  and 
has  to  be  determined  first  by  a step-by-step  procedure 

The  method  la  applicable  for  wings  with  arbitrary  subsonic  leading  or  trailing 
edges  and  for  wings  with  control  surfaces.  tn  tb t case  of  a sufe&onlc  leading  edge, 
the  upwash  field  of  boxes  cut  by  and  lying  in  front  of  the  leading  edge  are  presumed 
to  have  a square  root  singularity. 

The  formulation  of  present  M.B.B.  lifting  surface  theory  primarily  supplies 
the  velocity  potential  values  Sectional  or  total  loads  are  easily  derived  from 
these  values  by  applying  partial  integration  techniques 

If  the  calculation  of  pressure  coefficients  is  wanted,  this  can  be  achieved  by 
using  the  formula. - 


This  means  that  a numerically-given  function  has  to  be  dlffcrcn'  land  For 
this  reason,  the  function  must  be  smoothed  first  Caution  is  demanded  where  the 
potential  offers  a Hick,  l.e.  where  the  function  cannot  be  differentiated. 

6 . DISCUSSION  OF  RESULTS 

6 . I Comparison  of  Pressure  Distributions 

Theory  and  experiment  were  compared  for  four  xtreamvise  stations,  namely 
stations  1,  3.  5 and  6 (figures  3,  4.  5,  6>.  The  two  theories  agree  well  both  in 
their  real  and  imaginary  par!*  At  M • l 3,  the  B A C theory  predicts  pressures 
forward  of  the  nearl)  sonic  hinge  line,  this  being  due  to  the  smoothing  technique 
used  in  computing  pressure*  from  control  surface  nodes  Although  n all  cases, 
both  supersonic  theories  predict  pressure  maxima  in  the  proxlmit)  of  the  control 
surface  hinge  line,  the  maxima  obtained  from  wind  tunnel  teats  are,  in  gen*  ml. 
much  grester  This  is  probably  due  to  the  fact  that  neither  theory  accounts  for 
the  pressure  singularity  tn  the  case  of  a subsonic  hinge  line  <.n  all  cases,  the 
hinge  line  is  subsonic)  Along  station  6,  there  occur  vhe  greatest  dentations 
between  *est  and  theory,  with  a shif,  of  the  test  pressure  maximum  aft 

In  aoao  test  measurements  (e.£  stations  3 and  J> , figure  6).  there  appears  a 
secondary  maximum  of  pressure,  smaller  and  aft  of  the  hinge  line  peak 
This  Is  perhaps  evidence  of  shocks,  which  arc  not  considered  fcv  either 
t hco  r > . 

6 2 Coapar 1 son  of  Local Lift  and  Moment  Distributions 

I r.  figure  7.  loral  lift  and  moment  derivatives  plotted  virsus  span  are  sho»n 
The  parameters  are  M ? 1 l and  f ’ 145  Hz  corresponding  to  k 0 80 

The  supersonic  theories  agree  fatrl)  well,  but  the  comparison  between  thcor> 
and  experiment  does  not  show  satisfactory  correlation  in  all  cases,  especially  for 
the  out-of-phaae  derivatives  The  centre  of  pressure  is  further  outboard  tn  the 
tests  than  predicted  by  theory  and.  when  torsion  effects  are  present  in  an  aircraft 
flutter  analysis,  this  could  possibly  leaf,  to  larger  deviation* 

Figure  8 presents  sectional  lift  and  moment  distributions  for  M l 1 and 
f * 260  Hz  corresponding  to  k *■  1 40 

In  l n 1 » case,  there  is  good  agreement  between  the  two  theories  But  both 
predictions  overestimate  the  real  parts  of  lift  and  hinge  moment  at  the  inboard 
stations  and  underestimate  them  near  the  outboird  control  siue  edge  For  the 
imaginary  parts  of  lift  and  moment,  there  are  considerable  deviations  especially 
out  boa  rd 


The  data  for  U “ 1 3 and  f * 148  Hz.  corresponding  tr  k 0 70.  are  given  in 

figure  9 With  this  parameter  configuration,  we  have  deviating  predictions  for  the 
local  lift  distribution  whereas  M H B and  3 A C predictions  for  the  local  htnge 
moment  distribution  agree  well 


The  vxptrliv&t tl  cj  distribution  is  well  predicted  by  B.A.C.  theory  whereas  the 
distribution  c j**  is  well  predicted  by  M.B.B.  The  in-phase  hinge  moment  is  over- 
estimated  in  all  cases.  The  ou t -o f -phase  aoaent  coefficients  are  in  good  agreement 
with  test  valued 

Figure  10  presents  data  for  M - 1.3,  f * 260  Hz  and  k • 1.24.  As  in  the  p.evious 
graph,  the  agreement  of  the  experimental  spanwlse  lift  with  theory  depenas  on  the 
procedure  considered.  For  the  iasglnsry  coefficients  of  cj  and  cj..  we  -ran  observe 
alnor  differences  at  the  inboard  stations  and  greater  discrepancies  at  the  outboard 
station 

It  sees*  a feature  of  the  sectional  lift  that  theory  is  greater  than  test  near 
the  root  chord  and  vice-verse  at  the  tip.  If  the  planfora  is  not  sirro red  exactly 
which  is  probable,  due  to  the  presence  of  a boundary  layer  near  to  the  tunnel  wall, 
then  there  will  be  differences  in  results  between  test  sod  theory,  becoming  less 
significant  as  the  planfora  tip  is  spprosched,  since  the  theories  both  assume  s 
perfect  mirror  image.  On  station  1,  it  can  therefore  be  expected  that  theory  will 
produce  results  of  larger  magnitude  than  test,  and  this  is  shown  to  be  so  st  least 
for  real  parts.  Imaginary  parts  of  local  derivatives  are.  unfortunately,  much 
harder  to  predict  consistently.  This  is  fundamentally  due  to  their  being  the 
difference  between  two  relatively  large  opposing  effects. 

6.3  Comparison  of  Total  Lift  and  Total  Mince  Mone^*  Coefficients 


Figures  11.  12.  13  present  the  variation  of  the  total  lift  and  total  hinge 

moment  coefficient  with  reduewd  frequency.  Correlation  between  the  rcaults  from 
both  theoretical  methods  is  good.  This  is  not  cont radio tory  to  the  results  obtained 
from  comparisons  of  local  loads,  since  sectional  over-estimations  and  under- 
estimations  will  partially  cancel  each  other  when  Integrated  to  obtain  total 
coef  f iclents . 

Tbe  theoretical  values  of  lift  and  binge  moment  are  greater  in  magnitude  tha j 
those  of  experiment,  and  for  both  C ^ * and  Cjj*  , tbe  trends  with  k and  M are  correctly 
represented  Again,  it  is  found  that  the  Imaginary  parts  are  much  more  difficult  to 
predict  cons  * s t en 1 1 y . Ch“  is  particularly  critical  since  It  should  be  kept  positive, 
to  avoid  one  degree  of  freedom  flutter.  It  can  be  seen  from  figures  11-13  that  this 
is  so  for  all  cases  considered 

In  this  esse,  increasing  k is  favourable  and  increasing  M Is  unfavourable.  bo«b 
these  effects  being  predicted. 

7 CONCLUSIONS 

Both  theories  in  this  report  are  linearised  and  purely  supersonic,  whereas  the  tru' 
flow  is  mixed,  especially  at  tbe  lower  Mach  number.  The  Stark  theory  tends  to  prod*  e 
ntgbsr  pressure  magnitudes  than  Sa  ler-Allen.  and  a better  correlation  with  test  results 
This  1*  presumably  due  to  the  different  binge  line  treatments.  In  terms  of  overall 
forces,  however,  tbe  two  theories,  one  being  an  Integrated  dowpwasb  method,  tbe  other  in 
irtegrated  potential  method,  are  In  good  agreement 

In  the  context  of  aircraft  flutter  snaljsls.  therefore,  the  differences  between  these 
two  tbworeticsl  trcatmwnts  are  not  very  significant  In  generwl  . trends  demonstrated  oy 
experiment  are  followed  reasonably  by  both  theories,  but  it  is  of  greut  Importance  that 
both  theoretical  and  experimental  studies  are  of  very  high  accuracy,  so  that  control 
surface  flutter  statu*  can  be  established  with  confidence  st  low  supersrnlc  conditions 
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Fig.  5 Experimental  and  theoretical  pressures  at 
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